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ABSTRACT 

This paper describes the results of computational studies on a turbocharged four cylinder four 
stroke cycle direct injection diesel engine using Hydrogen as an alternative fuel to diesel. The 
simulation is done in the professional engine simulation software from AVL company of Austria 
named as BOOST. The modeling methodology involves the use of first law of thermodynamics 
and Navier-Stokes equations for engine cylinder and manifolds. The design parameters are fixed 
by engine geometry. A matrix was prepared for the operating variables to carry out the 
simulation. First the data was used as per diesel engine needs and results for engine performance 
were generated. The operation was revised with data for proposed hydrogen engine and its 
performance was studied. The software gave successful results in both the cases. It was observed 
that the power output and brake specific fuel consumption were improved with hydrogen as an 
alternative fuel due to its higher calorific value. It is proposed that hydrogen can be successfully 
used in diesel engine using pilot injection of diesel towards the end of the compression stroke or 
else by using a spark plug in diesel engine. It is proposed that only air will travel through 
compressor and intake manifold and liquefied hydrogen will be injected as per norms in multi 
point fuel injection systems using intake port for injection. 
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INTRODUCTION 

Computer simulation studies helps to predict the behavior of the engine in the diesel and 
hydrogen fuel modes. We prepare the models for Hydrogen and Diesel direct injection modes for 
the engine systems and feed the actual data corresponding to the design and operating conditions 
of the system. It helps to simulate the results without actually performing experiments. Thus a lot 
of money and time is saved. Moreover we can simulate and compute those results which are very 
difficult to be measured experimentally. Favorable computed results pave the way for further 
experimental investigations. 

The objectives are to investigate the feasibility of Hydrogen as future alternative fuel in diesel 
engines. It is also intended to choose hydrogen as fuel because it has no carbon atoms and least 
hydrogen atoms as compared to diesel and therefore will produce zero or minimal pollution. 
Table 1 at the end gives physico-chemical properties of hydrogen and diesel which help us to 
investigate the feasibility of using hydrogen as an alternative fuel to diesel. 
THEORETICAL BASIS 

The theoretical background including the basic equations for all elements used in the present 
model is summarized below to give a better understanding of the program. 
The Cylinder , High Pressure Cycle, Basic Equation. 

The calculation of the high pressure cycle of an internal combustion engine is based on the first 

law of thermodynamics: 

dirricu) pcdV dQp y> dQw hBB.dm.BB 
da da da da da (Eq.l) 



where 

d(mc.u) 
da 

pcdV 
da 

dQ F = 
da 



= change of the internal energy in the cylinder. 
= piston work, 
fuel heat input. 



sydQ^_ = -wall heat losses 
da 



A Monthly Double-Blind Peer Reviewed Refereed Open A ccess Interna tional e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., MJJiBBtffli^ as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Management, IT and Engineering 

http://www.ijmra.us 



331 



hBB.dlllBB tin n . n . 

= enthalpy flow due to blow-by 

da 

dniBB , , , „ 
= blow-by mass flow 

da 

The first law of thermodynamics for high pressure cycle states that the change of internal energy 
in the cylinder is equal to the sum of piston work, fuel heat input, wall heat losses and the 
enthalpy flow due to blow-by. 

Eq.l is valid for engines with internal and external mixture preparation. However the terms 
which take into account the change of gas composition due to combustion, are treated differently 
for internal and external mixture preparation. 
For internal mixture preparation it is assumed that 

The fuel added to the cylinder charge is immediately burnt. 

The combustion products mix instantaneously with the rest of cylinder charge and thus 
form a uniform mixture. 

As a consequence, the Air-Fuel ratio of the charge diminishes continuously from a high 
value at the start of combustion to the final value at the end of combustion. 
In order to solve this equation, models for the combustion process and the wall heat transfer, as 
well as the gas properties as a function of pressure, temperature, and gas composition are 
required. 

Together with the gas equation 

p c = _L .m c .R .T c (Eq.2) 

Establishing the relation between pressure, temperature and density, Eq. 2 for in-cylinder 
temperature can be solved using a Runge-Kutta method. Once the cylinder gas temperature is 
known, the cylinder gas pressure can be obtained from the gas equation. 
Combustion Model 

The following equation for the stoichiometric air requirement specifies how much air is required 
for a complete combustion of 1 kg fuel: 

L st = 137.85 .(— — +— — + — — ) [kg Air/kg Fuel] (Eq.3) 

12.01 4.032 32.06 32.0 
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For lean combustion, the total heat supplied during the cycle can be calculated from the amount 
of fuel in the cylinder and the lower heating value of the fuel. The lower heating value is a fuel 
property and can be calculated from the following formula: 

H u = 34835 . c +93870 . h +6280 . n +10465 . s -10800 . o -2440 . w [kj/kg] (Eq.4) 

In rich air fuel mixture combustion, the total heat supplied during the cycle is limited by the 
amount of air in the cylinder. The fuel is totally converted to combustion products even if the 
amount of air available is less than the amount of stoichiometric air. 
Heat Release Approach. 

The vibe function is used to approximate the actual heat release characteristics of an engine: 

^ = _^_ (m+1) y m e -a.y(m + l) (£q 5) 

da Aac 

dx = ^ (Eq.6) 

Q 

y = ot-- (Eq.7) 

Aac 

The integral of the vibe function gives the fraction of the fuel mass which was burned since the 
start of combustion: 

x= \(—.da) = l-e' dyim+1) (Eq.8) 

J da 

Gas Exchange Process , Basic Equation 

The equation for the simulation of the gas exchange process is also the first law of 
thermodynamics : 

d(m.u) p c .dV s-^ dQu ^ dm ^ dm 

— \ 1_ \ > (Eq.y) 

da da da da.hi da.he 

The variation of the mass in the cylinder can be calculated from the sum of the in-flowing and 

out-flowing masses: 

dm _ -yi dm y dm ^ ^ 

da da da 

Piston Motion 

Piston motion applies to both the high pressure cycle and the gas exchange process. 

For a standard crank train the piston motion as a function of the crank angle a can be written as: 
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s= (r+l).cos\|/-r.cos(\|/+a)-l. yi-{y.sin(^ + a)--} 2 (Eq.ll) 



\\i = arcsin( 



r + l 



) 



-(Eq.12) 



Heat Transfer 

The heat transfer to the walls of the combustion chamber, i.e. the cylinder head, the piston, and 
the cylinder liner, is calculated from: 

Q wi = Ai .a w . (T c -T wi ) (Eq.13) 

In the case of the liner wall temperature, the axial temperature variation between the piston TDC 
and BDC position is taken into account: 

i-e cx 



T, =T 



L,TDC 



X.C 



-(Eq.14) 



Tl.tdc 

c = ln { - — } 



-(Eq.15) 



Tl, BDC 

For the calculation of the heat transfer coefficient, the Woschni 1978 heat transfer model is 
used. 

Woschni Model 

The woschni model published in 1978 for the high pressure cycle is summarized as follows: 



1 on r»~0-2 0.8 rri -0.53 



L\.Cm + C2. .(p C - Pco) 



0.8 



-(Eq.16) 



PcuVc.l 

CI = 2.28+0.308.cu/cm 

C2 = 0.00324 for DI engines 

C2 = 0.00622 for IDI engines 
For the gas exchange process, the heat transfer coefficient is given by following equation: 

a w = l30.D-°- 2 .pc \Tc-°- 5 \(C3.c m ) 0S (Eq.17) 

C 3 = 6.18+0.417.c u /c m 
Fuel Injector 

The fuel injector model is based on the calculation algorithm of the flow restriction. This means 
that the air flow rate in the fuel injector depends on the pressure difference across the injector 
and is calculated using the specified flow coefficients. In addition, the amount of fuel specified is 
fed into the air flow. 
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In the case of carburetor model, the fuel flow is set to a specified percentage of the 
instantaneous mass flow. 

For the injector model, a measuring point must be specified at the location of the air flow meter. 
In this case the mean air flow at the air flow meter location during the last complete cycle is used 
to determine the amount of fuel. As is the case for continuous fuel injection, the fuelling rate is 
constant over crank angle. 

The fuel is added in gaseous form to the pipe flow. No evaporation is considered. 
Pipe Flow 

The one dimensional gas dynamics in a pipe are described by the continuity equation 

dp 8(pM) 1 dA m 18 > 

= r pu.—.— t (Eq.18) 

dt dx A ax 

the equation for the conservation of momentum 

= _ d( P M 2 + p)_ 2 l dA _ Fx ^ 9) 

dt dx A dx V 

and by the energy equation 

BE = _tfr<E + ,)] + I dA + 

dt dx A dx V 

The wall friction force can be determined from the wall friction factor Xf : 
FR Xf || 

= .pu\u\ (Eq.21) 

Using the Reynold's analogy, the wall heat flow in the pipe can be calculated from the friction 
force and the difference between wall temperature and gas temperature: 

Hi = -^L.p\u\.c P .(Tw-T) (Eq.22) 

During the course of numerical integration of the conservation laws defined in the Eq.20, Eq.2l 
and Eq.22, special attention should be focused on the control of the time step. In order to achieve 
a stable solution, the CFL criterion (stability criterion defined by Courant, Friedrichs and Lewy ) 
must be met: 

Ax 

At < — (Eq.23) 

u + a 

This means that a certain relation between the time step and the lengths of the cells must be met. 
The time step to cell size relation is determined at the beginning of the calculation on the basis of 
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the specified initial conditions in the pipes. However, the CFL criterion is checked every time 
step during the calculation. If the criterion is not met because of significantly changed flow 
conditions in the pipes, the time step is reduced automatically. 

An ENO scheme is used for the solution of the set of non-linear differential equations discussed 
above. The ENO scheme is based on a finite volume approach. This means that the solution at 
the end of the time step is obtained from the value at the beginning of the time step and from the 
fluxes over the cell borders. 
RESULTS AND DISCUSSION 
Effect of Speed on Power 

The Fig. 1 below shows the effect of speed on power. It is seen as the speed increases the power 
also increases due to more number of power cycles per unit time. Further it is seen that the power 
developed by hydrogen as fuel is much higher than diesel due to higher calorific value of 
hydrogen. 




4— — r- 

2000 4000 6000 

Speed, r.p.m. 
Fig.1 Effect of Speed on Power 



Effect of Speed on Brake Specific fuel Consumption. 

Fig. 2 below shows the effect of speed on BSFC ( brake specific fuel consumption ). It is seen 
that the operation of hydrogen engine is more economical as the BSFC of hydrogen fuel is less 
than that of diesel. It is due to higher calorific value of hydrogen. 
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Fig.2 Effect of Speed on BSFC 



Effect of Speed on Enthalpy Flow at entry to turbine. 

Fig. 3 below shows the effect of speed on enthalpy flow at entry to turbine in the selected 
turbocharged engine in cases of hydrogen and diesel as alternate fuels. Enthalpy flow at entry to 
turbine in case of hydrogen engine is higher due to better combustion and higher rate of heat 
released during combustion in cylinders. This will obviously result in larger enthalpy drop across 
the turbine in case of hydrogen fuel with the result more power can be recovered to run the 
compressor through turbine. 




Diesel 
Hydrogen 



2000 4000 6000 
Speed , r.p.m. 

Fig.3 Effect of Speed on Enthalpy 
Flow at entry to turbine 



A Monthly Double-Blind Peer Reviewed Refereed Open A ccess Interna tional e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., M.iJiBBtffi^ as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Management, IT and Engineering 

http://www.ijmra.us 



Effect Of Air-Fuel Ratio On Power 

The Fig. 4 shows the effect of Air-Fuel ratio on power. It is seen as the the air-fuel ratio 
increases the power output decreases as the mixture becomes leaner. Further it is seen that the 
power developed by hydrogen as fuel is much higher than diesel due to higher calorific value of 
hydrogen. Highest power is developed at an air-fuel ratio of about 20 in both cases. This is true 
for diesel , as diesel engines are operated on this air-fuel ratio for maximum power generation. 
The air-fuel ratio for diesel engines cannot be decreased further as it produces lot of smoke due 
to heterogeneous combustion. It is seen that hydrogen engines can also be operated on an air-fuel 
ratio of about 20 for maximum power generation. 



250 




20 40 60 

Air-Fuel Ratio, (-) 

Fig.4 Effect of Air-Fuel Ratio on Power 



Effect Of Air-Fuel Ratio On Brake Specific Fuel Consumption.(BSFC) 

Fig. 5 shows the effect of air- fuel ratio on BSFC ( brake specific fuel consumption ). It is seen 
that the operation of hydrogen engine is more economical as the BSFC of hydrogen fuel is less 
than that of diesel. It is due to higher calorific value of hydrogen. Further hydrogen engine 
remains economical over a wide range of air- fuel ratio. This feature will make hydrogen a more 
suitable fuel for automotives and power generators over a wide range of speed and load. 
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Effect Of Air-Fuel Ratio On Enthalpy Flow At Entry To Turbine 

Fig. 6 shows the effect of air-fuel ratio on enthalpy flow at entry to turbine in the selected 
turbocharged engine in cases of hydrogen and diesel as alternate fuels. Enthalpy flow at entry to 
turbine in case of hydrogen engine is higher due to better combustion and higher rate of heat 
released during combustion in cylinders. This will obviously result in larger enthalpy drop across 
the turbine in case of hydrogen fuel with the result more power can be recovered to run the 
compressor through turbine. This will also help to increase the power and thermal efficiency of 
hydrogen engine through turbo-compounding. 
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Fig.6 Effect of Air-Fuel Ratio on 
Enthalpy Flow at entry to turbine 



CONCLUSIONS 

a. Hydrogen can safely be used in diesel engines as an alternative fuel to conventional diesel. 

b. Pollution formation from engines using Hydrogen as fuel will be minimal due to no carbon 
atoms in hydrogen fuel. 

c. Power is increased by using Hydrogen as fuel in diesel engines. 

d. The Hydrogen engine operation is economical as the brake specific fuel consumption with 
Hydrogen is less as compared to diesel. 
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APPENDIX-A 
NOMENCLATURE 



a = speed of sound 

A = pipe cross-section 

Aeff = effective flow area 

A; = surface area (cylinder head, piston, liner) 

AFcp = air fuel ratio of combustion products 

Ageo = geometrical flow area 

c = mass fraction of carbon in the fuel 

c v = specific heat at constant volume 

c p = specific heat at constant pressure 

CI = 2.28+0.308.cu/cm 

C2 = 0.00324 for DI engines 

C2 = 0.00622 for IDI engines 

Cm = mean piston speed 

Cu = circumferential velocity 

c u = circumferential velocity 

D = cylinder bore 

D = pipe diameter 

dm; = mass element flowing into the cylinder 

dm e = mass element flowing out of the cylinder 

d v i = inner valve seat diameter (reference diameter) 

^ mBB = blow-by mass flow 
da 

e = piston pin offset 

1 , 

E = energy content of the gas (=p. cv T + — .p.u ) 

f = fraction of evaporation heat from the cylinder charge 

Fr = wall friction force 

h = mass fraction of hydrogen in the fuel 

h B B = enthalpy of blow-by 
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h; = enthalpy of in-flowing mass 

h e = enthalpy of the mass leaving the cylinder 

H u = lower heating value 

k = ratio of specific heats 

1 = con-rod length 

m = shape factor 

m = mass flow rate 

m c = mass in the cylinder 

m ev = evaporating fuel 

m p i = mass in the plenum 

n = mass fraction of nitrogen in the fuel 

o = mass fraction of oxygen in the fuel 

p = static pressure 

Poi = upstream stagnation pressure 

Pc,o = cylinder pressure of the motored engine[bar] 

Pc,l = pressure in the cylinder at rVC[bar] 

p p i = pressure in the plenum 

p c = cylinder pressure 

P2 = downstream static pressure 

q ev = evaporation heat of the fuel 

q w = wall heat flow 

Q = total fuel heat input 

Qf = fuel energy 

Q w i = wall heat flow (cylinder head, piston, liner) 

r = crank radius 

Ro = gas constant 

s = piston distance from TDC 

t = time 

T = temperature 

Tc,l = temperature in the cylinder at intake valve closing (IVC) 

T c = gas temperature in the cylinder 
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T w ; = wall temperature ( cylinder head, piston, liner) 

T L = liner temperature 

T l,tdc = liner temperature at TDC position 

T l,bdc = liner temperature at BDC position 

T w = pipe wall temperature 

Tqi = upstream stagnation temperature 

u = specific internal energy 

u = flow velocity 

V = cylinder volume 

V = cell volume (A.dx) 

VD = displacement per cylinder 

w = mass fraction of water in the fuel 

x = relative stroke (actual piston position related to full stroke) 

x = coordinate along the pipe axis 

a = crank angle 

a = start of combustion 

Aa c = combustion duration 

a w = heat transfer coefficient 

p = density 

urj ■ 1= flow coefficient of the port 
v|/ = crank angle between vertical crank position and piston TDC position 

Xf H = wall friction coefficient 

At = time step 

Ax = cell length 
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APPENDIX-B 



Turbo charged Diesel Engine Specifications 


Bore 


90 mm 


Stroke 


98 mm 


Compression Ratio 


19 


Number of Cylinders 


4 



APPENDIX-C 

Table 1: Physico-Chemical Properties of Hydrogen and Diesel. [2] 



Fuel Property 


Hydrogen 


Diesel 


Formula 


H2 


C8 to c25 


Molecular weight 


2.02 


200 approx. 


Composition, weight% 






Carbon 





84-87 


Hydrogen 
Oxygen 


100 



13-16 



Freezing point 


-275 


-40 to -1 


Lower heating value, Kj/g 


119.93 


42.5 


Flammability limits,vol% 






Lower 






Higher 


49 


1.0 

6.0 




75 


Stoichiometric air-fuel ratio, 


34.3 


14.7 


weight 






Cetane number 




40-55 
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